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Abstract. We investigate the orbital stability and instability of standing 
waves for two classes of Klein-Gordon equations in the semi-classical regime. 



1. Introduction and results 

The nonlinear Klein-Gordon equation 

(1.1) £ 2 u tt -e 2 Aw + mu - \u\ p ~ 1 u = (t, x) G R + x R N , 

where e,m > 0, p > 1 for N = 1,2 and 1 < p < (N + 2)/(iV - 2) for N > 3, is one 
of the simplest nonlinear partial differential equations invariant for the Poincare 
group. We are interested in the study of the nonlinear Klein Gordon equation in 
presence of a potential depending on the space variable. Two different choices are 
viable. We can simply add a potential term W{x)u to equation (1.1). This case has 
been studied, for the linear wave equation, for example, by Beals and Strauss in [3]. 
Otherwise, if we want to fully preserve the invariance for the Poincare group of 

(1.1) , we have to change the temporal derivative e 2 dtt with a covariant derivative, 
depending on the potential D 2 t , where D t = edt + iV(x). This second approach 
is classical, in quantum electrodynamics, when considering electromagnetic waves. 
The first approach leads us to consider the equation 

(1.2) £ 2 u tt -E 2 Au + mu-Wu-\u\ p - 1 u = 0, in R N , 
while the second one to investigate the equation 

(1.3) s 2 u tt + 2ieVu t - s 2 Au + mu-V 2 u- \u\ p ~ 1 u = 0, in R N . 

In this paper, we treat simultaneously the two previous Klein-Gordon equations by 
studying 

(1.4) e 2 u tt + 2ieVu t - e 2 Au + mu- Wu~ = 0, in R N , 

where u : R N x M -> C, e > and V, W are real valued potential functions. 
Equation (1.4) formally yields to (1.3) for the choice W = V 2 as well as to (1.2) 
when V — 0. We shall study the stability of standing waves of this equation in the 
semiclassical regime e — > 0. It admits standing waves, namely solutions of the form 
u(x,t) = e lu}t / e {p u (x/e), where wel and ip u satisfies 

(1.5) - A<p u + (m -uj 2 - 2loV(e V ) - W{ey))ip u - |^| P "V^ = 0, in R N . 
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To ensure existence of solutions to (1.5) for e close to 0, we assume the following. 
Let V and W be C 2 . For the function 

Z(y) = m -uj 2 - 2ujV(y) - W(y), y £R N 

there exists xo G R w such that 

(1.6) VZ(xo) = 0, S7 2 Z(xo) is non-degenerate. 
Furthermore, we assume that 

(1.7) Z(x ) =m-uj 2 - 2ujV(x ) - W{x ) > 0. 

Under these hypotheses, it is well- know (sec e.g. [1] or [2, Section 8.2]) that 
when e is close to the equation (1.5) admits a family of positive, exponentially 
decaying, solutions ip u C _ff 1 (]R JV ) (hiding the dependance upon e). More precisely, 
there exists £ e € and ip^ G H 1 (R N ) such that tp u (-) = ip u (- - f e ) + 0{e 2 ) in 
H 1 ^) as e -> 0, where £ e = o(e) and 

-A^ + Z(x Q )^ = \il>u\v- l i> u , in R N . 

Wc are interested in the (orbital) stability or instability of standing waves when 
e goes to 0. 

A standing wave of (1.4) is said to be (orbitally) stable if any solution of (1.4) 
starting close to the standing wave remains close for all time, up to the invariances 
of the equation. More precisely, we say that e~cp u (^ j s stable if for all i] > 
there exists 6 > such that for all uq € i/ 1 (R Ar ) verifying ||«o — (Pu>\\h>-(r n ) < 8 
the solution u(x,t) of (1.4) with initial data uo satisfies 

sup inf ||u - e 4 V w ||.ffi( R jv) < V- 
ten eeR 

Since the pioneering works [4,6-8,24,25], the study of orbital stability for standing 
waves of dispersive PDE has attracted a lot of attention. Among many others, one 
can refer to [10,11,13]; see also the books and surveys [5,12,22,23] and the references 
therein. Relatively few works [9, 14, 16] are concerned with stability at the semi- 
classical limit for Schrbdinger type equations. For Klein-Gordon equations, after 
the ground works [20,21], there has been a recent interest for stability by blow- 
up [15,17-19]. 

To study stability, we first rewrite (1.4) in Hamiltonian form 

(1.8) = JE'(U), 

where U=(fj,J= ^ , and 

E (U) = ^\\v - iVu\\ 2 L2{wN) + -\\Vu\\ 2 L2{mN) + ml\\u\\ 2 L2(RN) 



-\ \ W \u\ 2 dx-^-\\u\\l + p 
2 Jrn p + 1 

It is easy to see that if u solves (1.4) and v is defined by v := eu t + iVu, then 
U = I ) solves (1.8). The charge Q(U) is defined by 



Q(U) = Sj / uvdx. 
Jr n 
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In particular, for a standing wave u = e llJjt ' e (p u (x/e), the charge is given by 
(1.9) Q{^):=Q(U)=e N L\\^\\ 2 LHm) + f V(ey)\<p u 

According to the theory developed in [7,8], a standing wave e^ipa, (|) is stable 
if two conditions are satisfied. 

d 

(i) The Slope Condition: -~-Q(tpJ) < 0. 

U(jJ 

(ii) The Spectral Condition: L £ := —A + Z(ey) — p|v 3 oj|' p_1 has exactly one 
negative eigenvalue and is non degenerate. 

On the other hand, denote by n(L e ) the number of negative eigenvalues of L e and 

d d 
set p(w) = if — Q(<Pu) > 0, p(uj) = 1 if —Q(tp u ) < 0. Then the standing wave 

is unstable if n(L e ) — p(u>) is odd. 
Then we have the following 

Theorem 1.1. Assume that conditions (1.6)-(1.7) hold. Then, we have the follow- 
ing facts. 

(1) If p < 1 + i/N, then the Slope Condition t£jQ(<£> w ) < is fulfilled if 

Z(xo) < (uj + V(xoj)' 2 ( N] (non-critical case) 

\v — 1 / 



or if 



Z(x ) = (oj + V(x )) 2 (^-n), 

I or iwim^ (critical case). 

AZ(0) - AF(0) (l + 2( %p) )) < °i 

(2) If p > 1 + i/N, then we always have -^jQif^) > 0. 

(3) We have the equality n(L s ) = n(V 2 Z (xq)) + 1, where n(V 2 Z(xo)) is the 
number of negative eigenvalues o/V 2 Z(xo). 

In particular, the standing waves e lLJt Lp u are stable if xq is non-degenerate local 
minimum of Z , p < 1 + A/N and 

Z(x ) <(lu + V(x )) 2 (-^ - N) . 

Note that, conversely to what was happening in the case of Schrddinger equations 
studied in [14], the values of the potentials V and W at xq comes into play for the 
Slope Condition even in the noncritical case. Note also that only the local behavior 
of Z around xq influences the stability or instability. 

Notations : Most of the objects we consider will depend both on e and ui. We 
will emphasize the most important parameter by indicating it as a subscript, the 
dependence in the other parameter being understood. 

2. Proof of Theorem 1.1 

In this section, we prove Theorem 1.1. We start be focusing on the Slope Con- 
dition and then we study the Spectral Condition. For the sake of simplicity in 
notations and without loss of generality, in the rest of this section we assume that 
xq = 0. 

2.1. The Slope Condition. We start with the noncritical case. 
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2.1.1. Noncritical case. We assume that 

Z(0)^( W + F(0)) 2 (^i-iv). 

We first rewrite Q(<Pu) by expanding V{ey) and using the exponential decay of 
Q(<p u ) =e N (u; + V(0)) ||^||^ (RW) + 0(e N+1 ). 

Therefore, since 

(2.1) £-Qfr u ) = e N yJl 2{MN} + e N (u> + V(0))-^-yJl 2{mN) + 0(s N+1 ), 



to evaluate the sign of the map ui i-> ■M-:Q( t Pu>) one should compute the quantity 



(2-2) ^ll^|il,2 (R iv ) = 2 y ^ iC<A^, 

where := (x). 

We remark that differentiation of (1.5) with respect to w easily yields 

(2.3) L B R u =2((j + V(ey))( Pa . 

If we now introduce the rescaling (f u (x) = A^ 1 <p x (\f\x), it follows that ip x satisfies 



(2.4) — Aif \ + X^ 1 Z i —j= j tp\ — |^rVA = 0, 



1 l,„. IP-1 n :„ TD>iV 

Now, differentiating equation (2.4) with respect to A and denoting T x = ^§ x L \ x _ l 
yields 

(2.5) L £ T X - Z(ey)tp u - X -ey ■ VZ(ey)^ u = 0. 
Since is a critical point of Z , a Taylor expansion gives 

(2.6) Z(ey) = Z(Q) + 0(e 2 \y\ 2 ), 

(2.7) \eyVZ(ey) = 0{e 2 \y\ 2 ). 
Then, from (2.5), as e — > we have 

(2.8) L £ T x = Z(0)tp u + O(e 2 \y 2 \tp u ), in R N . 
Then, in turn, taking into account identity (2.3) we get 

Z(0) [ R u ip u = [ R^L e T\ + 0(e 2 ) 
= [ L^Tx + Oie 2 ) 

(2.9) = / 2(w + V(ey))tp u T x + 0(s 2 ) 



2(u J + V(0)) / <p u T x + 0(e) 

(u + V(0))^\\v x \\1 2(rN)ix=i +O(e) 

(w + V(0))(y - ^j) ||^|| 2 L2(KN) + 0(e). 
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In conclusion, by combining (2.1), (2.2) and (2.8), we have 
duj 



Then, taking into account the fact that Z(0) > and that converges to ip u in 

du> ^ 

4 



L 2 (R N ) as e — > 0, the sign of ^r-Q(ip u ) is the sign of 



z(o)-(c + ^(o)) 2 (— t -jv). 

2.1.2. Critical case. We assume now that 

(2.10) z(0) = (u + V(0)) 2 (^-N). 

In the critical case, the term of order e N in the expansion of -J^Q(cp u ) vanishes and 
we need to calculate the expansion at a higher order. We first refine (2.6)-(2.7). 

Z(ey) = Z(0) + y V 2 Z(0)(2/, y) + 0(e 3 \y\ 3 ) 
l -ey ■ VZ(ey) = E - V 2 Z(0)(y, y) + 0(s 3 \y\ 3 ). 

Then (2.5) gives 

L e T x = Z{Q)ip u + £ 2 V 2 Z(0)(y, y)<p u + 0(e 3 \y 3 \)^- 

Now, we have 

(2.11) Z(0) [ R u tp u = [ R^L £ T x -e 2 f V 2 Z(0)(y, y)R uVu + 0{e 3 ). 
From (2.3), we obtain 

(2.12) f R u L e T\ = [ L e R u T x = [ 2(cj + V(ey))ip u T\. 

JR N it" Jl» 

Expanding the potential V we get 

(2.13) / 2V(ey)(p u T x = f 2V(0)<p u T x + 2e f y V7(0)^T A 

Jm. N Jm. N JR N 

+ e 2 [ V 2 y(0)(y,2/)^T A + O( £ 3 ). 



/r n 

Note that since = i/) u (- — £ £ ) + C(e 2 ) and £ £ = o(e), we have 

(2.14) 2e f y W(0)^T A =2e f y VV(0)^T x + o{e 2 ) = o(e 2 ) 

where the last cancellation comes from the fact that ip u is radial. Coming back to 
(2.12) and as in (2.9), we have 

(2.15) / J R w L £ T A = ( W + F(0))(^-^-)||^|| 2 2(MN) 

+ e 2 f V 2 V(0)(y,y)<p u T x + o(s 2 ). 



It remains to compute the integrals involving the Hessians in (2.11) and (2.15). 
Since our problem is invariant by an orthonormal transformation, we can assume 
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without loss of generality that V 2 F(0) = diag(6i, . . . , 6jv)- Hence \7 2 V(0)(y, y) 



2j=i bjyj. Recall also that T\ can be computed explicitly to have 



Therefore, 



T\ = ^— (flu - \y ■ V(p u 

p — 1 2 



dyk ' 



We have after integration by parts 

N n N 



2 2 



2 E/ yjV^^r = -Hl (y* + 25 jk y])vi = -(N + 2) [ 
£r{Js. N dy k t^Jui" Jm 1 

Thus 

/ ^V(0)(y,y)^T x = f2 [ bj y^ u T x = - (J- - *±*) f) 6,/ 

Recall the following expansion in e for and yj^. 

<Pu = 4>u + o(e), i? w = ^^+o(e). 

aw 

Therefore, since i/'w is radial, 

/ 2/|^+o(e) = ^||M^||i2 (R « ) +o(e), 

JR N ly 

and so 



2 2 



(2.16) / VY(0)( H )^T A = 

1 7V + 2\ 1 



p- 1 4 y TV 



|y|^||i 2(RN) Al/(0)+ O (e). 



Similarly, we have 

(2.17) / V 2 Z(0)(y,y)i?^ = 



HllM^I!i2 (R «)Az(0) + o ( £ ). 



j>- 1 4 y n' 

Summarizing, using successively (2.11), (2.15), (2.16), (2.17) and (2.10) we have 
obtained 



(2.18) 1^ = - 2(CJ+ V (Q)) II^||| 2( ^ ) 

2 (AZ(0) - AV(0)) (I N + 2 



) IllylV'ullia^) +o(e 2 ). 



7VZ(0) VP" 1 

Now, we compute ' . First, recall that, coming back to the definition (1.9) 

of Q, we have 

-N d Q(<P«>l 

duj 



ON THE STABILITY OF KLEIN-GORDON EQUATIONS 



7 



As in (2.13), (2.14), and (2.16) we can expand in e and get 
V{ey)R u <p u = 27(0) / R u <p u 



Therefore, 



Using (2.18), we finally get 



U - 1 4 J N 



2 



i Z( 0)-AV(0)(l + M?»)) + o«. 

2.2. The Spectral Condition. We define the operator L Q := -A + Z(0) -pipZ~ l . 
It is well known (see e.g. [2]) that the spectrum of Lq consists of one negative 
eigenvalue, a iV-dimensional kernel (generated by for j = 1, . . . ,N) and the 
rest of the spectrum is bounded away from 0. When e is close to 0, the spectrum of 
L e will be close to the spectrum of Lq. In particular, the eigenvalue, of multiplicity 
N, will transform into N possibly different eigenvalues close to but shifted cither 
to the positive or to the negative side of the real axis, depending on the sign of the 
eigenvalues of the Hessian of Z at 0. More precisely, the following proposition was 
proved in [14] (see [9] for a detailed justification). 

Proposition 2.1. The spectrum of L £ consists of positive spectrum away from 
and a set of N + 1 simple eigenvalues {Aq, Ai, . . . , Xn} such that 



A < Ai < ■ ■ • < A 



N- 



As e — > 0, we have Xq < and the following asymptotic expansion holds for the 
other eigenvalues: 



Xj = Cj e 2 + o{e 2 ), j = l,...,N, 



where Cj = ^ ^„ |' 2 2( * Ml aj and {ai, . . . , on} are the eigenvalues of the Hessian 
matrix V 2 Z(0). 

Therefore, (3) in Theorem 1.1 is a direct consequence of Proposition 2.1. In 
particular, the spectral condition for stability will be satisfied if and only if is a 
non-degenerate local minimum of Z. 
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